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I. Executive Summary 

This report presents an analysis of the sound spectra generated by a trailing aircraft 
vortex during its rolling-up process. The study demonstrates that a rolling-up vortex 
could produce low frequency (less than 100 Hz) sound with very high intensity (60 dB 
above threshold of human hearing) at a distance of 200 ft from the vortex core. The 
spectrum then drops off rapidly thereafter. 

A rigorous analytical approach has been adopted in this report to derive the spectrum 
of vortex sound. First, the sound pressure was solved from an alternative treatment of the 
Lighthill’s acoustic analogy approach [1], After the application of Green’s function for free 
space, a tensor analysis was applied to permit the removal of the source term singularity of 
the wave equation in the far field. Consequently, the sound pressure is expressed in terms 
of the retarded time that indicates the time history and spacial distribution of the sound 
source. The Fourier transformation is then applied to the sound pressure to compute its 
spectrum. As a result, the Fourier transformation greatly simplifies the expression of the 
vortex sound pressure involving the retarded time, so that the numerical computation is 
applicable with ease for axisymmetric line vortices during the rolling-up process. 

The vortex model assumes that the vortex circulation is proportional to the time and 
the core radius is a constant. In addition, the velocity profile is assumed to be self-similar 
along the aircraft flight path, so that a benchmark vortex velocity profile can be devised 
to obtain a closed form solution, which is then used to validate the numerical calculations 
for other more realistic vortex profiles for which no closed form solutions are available. 

The study suggests that acoustic sensors operating at low frequency band could be 
profitably deployed for detecting the vortex sound during the rolling-up process. 
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II. Introduction 

The study of aircraft vortices has been a subject of both academic and practical interests 
for decades [2] , [3] . Some of the comparatively matured techniques for characterizing and 
tracking wake vortices from full scale aircraft have been flow visualization (via wing-tip 
smoke injection, smoke screen or natural condensation), in-flight probes on penetrating 
aircraft, instrumented fly-by towers, propeller and sonic anomometer arrays, acoustic 
radars, LDV/CW-Lidar and pulsed-Lidars [4]— [15] . Exploratory efforts have also been 
made in the areas of, for examples, RASS, Radar, pressure transducer array, infrared 
sensing, microwave radiometry and passive acoustics [3], [4], [16]— [20] . The present report 
examines the acoustic emission of aircraft vortices, in response to the renewed interest in 
studying the passive wake acoustic detection scheme. 

A number of held observations have indicated that vortices emit audible sound when in 
the near proximity of the ground. Frequencies from 1 to 5 KHz have been reported [4] , and 
various descriptions of the wake sound have been described from “whooshing” , “whine” , 
“crackling” to “roaring” [4], [7]. During early days of wake vortex sensor research, the 
aforementioned acoustic emission phenomenon was not deemed reliable enough as the 
basis of a sensor for tracking wakes near the ground. Consequently, no effort had been 
devoted over the years to systematically study the acoustic emission from aircraft wakes. 
More recently, unpublished measurements of exploratory nature from NOAA [21] indi- 
cated that vortices also emit infrasound. The general area of aircraft wake acoustics 
however, unlike other branches of aeroacoustics involving jets, airframe and blade-vortex 
interaction, has not received the comparable level of attention over the course of noise 
research. Consequently, fundamental issues such as the spectral characterization, its as- 
sociated uniqueness and consistency, as well as over what frequency band could wake 
vortices be tracked effectively and practically, still remain to be addressed. In response 
to the renewed interest in tracking and characterizing wakes with a passive acoustics ap- 
proach, the present report represents a first step in analytically examining the spectral 
characteristics of aircraft vortices during its initial generation phase at altitude. 

The analytical study of sound generation aerodynamically has been a continuing re- 
search area since the acoustic analogy approach pioneered by Lighthill [1], Representa- 
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tive studies that have since been influenced by Lighthill, examined the acoustic signature 
stemming from stretching or spinning vortex rings [22], and by ring vortices interacting 
with structure [23] or with other vortex rings [24], In these studies, the vortex velocity 
was treated as a source and the vortex pressure was evaluated by integrating the product 
of the time-retarded Green’s function and the source term over the entire region of vor- 
tex velocity distribution [23], [25]. The time-retarded Green’s function was expanded in 
terms of its retarded time in the far field, so that the vortex pressure was expressed in the 
form of a Taylor series and evaluated by truncating higher order terms. In studying the 
sound generated by aircraft wakes, the present report adopted a tensor analysis instead, 
to simplify the “source” term in the acoustic analogy approach. The spectrum of vortex 
pressure is expressed as the trace of the integration of a matrix. The study then examines 
the acoustic spectrum of a straight line vortex at altitude developing over a microphone 
array. A closed form solution of the spectrum of the vortex pressure field is obtained for 
a benchmark vortex profile, and is used for validation of the numerical method for other 
profiles where closed form solutions are not available. The development of the analysis is 
presented as follows. 

III. Formulation of the Governing Equation 

A. Vortex Sound Wave Equation 

The sound pressure, p, of an inviscid flow field is governed by the following set of 
equations [23] 


+ V • (pu) = 0 

(1) 

d 

— (pu) + V • (pw w) + Vp = 0 

(2) 

dp = c 2 0 dp 

(3) 


where c Q is the sound speed, p is density, p is pressure, and u is velocity field. Following the 
arguments of Kambe [23] for an incompressible flow, i.e., p = p Q the equations of motion 
are reduced to a wave equation involving pressure, classically known as the Lighthill’s 
acoustic analog equation [1]: 

1 0 ^ 

V ' 2 P - # opp = -p ° v ' ^ 


( 4 ) 
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The right hand side of (4) may then be interpreted as a “source” term for generating the 
pressure p. Far away from a source where the velocity u is zero or a constant, Eq. (4) is 
then satisfied by a spherical wave solution. 

B. Solution of Vortex Sound Pressure 

The solution of the vortex sound pressure from the wave equation (4) can be expressed 
as an integration involving the Green’s function as [26] 

p(r,t) = —p 0 JJ g(f — f';t — t'){v' • V' • [u{f\ t')u{f\ t')]^df'dt' (5) 

where g is the time-retarded Green’s function in free space 

0 if t < t' 


g(r — r';t — t') = 


—5 [|r — r'\ — c Q (t — t')\ if t > t' 


( 6 ) 


47t | r — r' 

with 5 representing the delta function. By applying the identities described in Eqs. (49) 
and (51) in Appendix A, and using the divergence theorem where the surface integral 
terms are ignored, one obtains 


p(r, t ) = — p 0 trace / / (u(r', t')u(r', t')) ■ V'V'g (r — r'; t — t') dr'dt' 


The expression with the index notation is then obtained as 

P(r, t) = ~Po JJ t')ui(f', t ')) (r ~ r'] t - t') dr'dt ' 

^ J 


Note that the relation of 


d d d d 


dx\ dx'j dxi dxj 
holds for the Green’s function, then Eq. (8) becomes 

d d 


p(r,t) = -p 0 — — JJ (ui(r', t')uj(r', t')) g (r — r']t — t 1 ) dr'dt' 


dxi dxj 

Placing the Green’s function into Eq. (10) yields 

d d 


p(r , t) = c 0 p 0 


dxi dxj 


{uiir' ,t')uj(r' ,t')) 


6 r — r'\ 

- c Q (t - t')' 

47 r 

\r - r'\ 



dr'dt ' 


( 7 ) 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


Furthermore, by performing the integration over t', one obtains 


P(r,t) = p c 


d d r 1 ( , \r — r'|\ ( , |r — r'\ 


dxi dx; J 477 | r — r 


7i Ui r', t — 


Uj r', t — 


Co 


dr' (12) 



Applying the property 


(13) 


a_ 

dxi 


Xi — x\ d 


\r — r'\ c 0 dt 

in the far held where |r — f'\ « r and ( Xi — x[) ~ Xi, Eq. (12) is then approximated as 

p 0 XiXj d 2 


P(r,t) 


47t rc 2 r 2 dt 2 


Ui r , t 


r — r 


Uj r',t 


I /y* I 


dr 


(14) 


Note that Eq. (14) is identical to the expression provided by Powell in Ref. [27]. 


C. Remarks on Retarded Time 


The retarded time, t — \f — f'\/c 0) that appears in Eq. (14) where c Q is the speed of 
sound, is introduced by the retarded-time Green’s function (6). It describes the time 
history of the vortex sound sources. In the far-held, where |f — f'\ « r — (r/r) • f r , the 
retarded time is approximated as 


t — 




(15) 


Note that the last term on the right-hand side of the above equation is the function of the 
orientation for the source position (?) with respect to the direction of observation (f/r), 
and it does not depend on the distance to the observer (since f/r is a unit vector). Thus 
this term should ideally be kept when carrying out the integration over r' as performed 
in Eq. (14), which formed the basis of subsequent development. 


IV. Formulation of Vortex Sound Spectrum 

In obtaining the spectrum of a vortex sound pressure, the Fourier transformation can be 
applied directly to Eq. (5) without performing complicated mathematical manipulations 
from Eq. (7) to Eq. (14). Denoting the term in the integrand as 


f(f' 1 t')=V-[V-u(f' 1 t')u(f\t')} 


(16) 


and utilizing the selection property of the delta (6) function, a more compact form for 
Eq. (5) is then obtained, 




r, t — 


\r — r 


df' 


( 17 ) 
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Note that Eq. (17) represents a general solution of vortex pressure field governed by 
Eq. (4) for an arbitrary profile of velocity distribution. Theoretically, knowing the veloc- 
ity profile u, the pressure can be evaluated directly, although the calculation involving 
double gradients on the velocity field may be sensitive to errors introduced by numerical 
procedure. Next, to obtain the spectrum of the vortex sound, Fourier transformation is 
applied directly on the integral representation of the pressure, and then an alternative 
form is obtained by operating the double gradients on the scalar Green’s function, instead 
of the vortex velocity field. 


A. Fourier Spectrum of Vortex Sound Pressure 

The spectral character of vortex sound is obtained by the Fourier transformation of 


P(r,u>) = 


p ° 




(18) 


VK\r-r'\ \ c 0 

where u is angular frequency. By writing the retarded time as t' = t— \ f — f'\/c 0 , Eq. (18) 
becomes 

r P Mr-r'\/co 

(19) 


■7 / , d f ■ r" f I / /■» 

P(f,uj) = Po I — —F{r',uj)df' 


4-7T | r — r 

where F{r',u ) is the Fourier transform of f(f',t'). In addition, defining the tensor 

A(r',u) = f [u(r', t')u(f', t')] e juJt 'dt' 


and writing the scalar Green’s function as 


(f>(r, r 1 ) = 


e jL0\r-r’\/c o 

47T If — f | 


( 20 ) 


( 21 ) 


Eq. (19) is simplified as 


P(r,u>) = Po trace / (A(r',u) ■ V'V'((i(r,r') dr 


( 22 ) 


The details of the derivation are found in Appendix A. The trace operator is the sum- 
mation of diagonal elements in the matrix representation of a tensor. Given the vortex 
velocity profile u, the tensor A can be calculated from Eq. (20). In addition, V'V'i j> can 
be computed using the definition of (j) in Eq. (21). Finally, the spectrum of vortex sound 
pressure expressed in Eq. (22) can be evaluated quantitatively. Note that the evaluation 
of the spectrum involves integrations over time t' (20) as well as local space r 1 (22). 
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B. Far-Field Representation of Vortex Sound Spectrum 

For practical reasons, far-held approximation of the pressure spectrum for vortex sound 
is of particular interest. In the far held, where r 3> r' and |r — f'\ ~ r — r ■ r’/r, the term 
V'V'0 in the integrand of Eq. (22) becomes, by keeping the lowest-order term with respect 
to 1/ | f — f' |, 

u 2 e jui\r-r'\/c 0 

VVf « -ff 9 .. (23) 

c- 7 t r 6 

Therefore the spectrum in the far held can be expressed approximately as 

P(r,uj ) ~ —trace { %, — ° e^ r ^ c ° [ A(r',u > ) • ^ e _J ^ r ' r df'l (24) 

( Cp7rr J r 2 J 

where cu 2 is due to the double-gradient operation of the scalar Green’s function. Note 
that the spectrum strength is inversely proportional to the distance r by knowing that 
(r — f') (f — f')/r 2 approaches to a unit tensor for large r. Also note that, at u = 0, the 
pressure spectrum is zero, i.e., P(f,co = 0) = 0, which confirms that vortex sound does 
not contain a DC component. 

Given the velocity prohle, the general form of spectrum for vortex sound in the far held, 
Eq. (24), can be evaluated numerically. Axisymmetry of the vortex velocity is assumed 
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for the analysis, whose nomenclatures is illustrated in Fig. (1). The axisymmetric line 
vortex field in space and time is also assumed to have the following general form, 

u (r, if) = (j) [S' (z! - Li) - S (V - L 2 )\ (p, t') (25) 

where (j) 1 is a unit vector in the circumferential direction, and S is the step function defined 



and <£>( p',t ') is the scalar velocity field. 


z <0 
z > 0 


(26) 


V. Rolling-Up Vortices with Constant Longitudinal Distribution 

The vortex during the rolling- up process is assumed to be longitudinally constant, so 
that the variables p' and t' are separable and the velocity is re-written as 

u (f\ 0 = $ [S (/ - t,) - 5 y - i 2 )] 4 (p') /((') (27) 

where f(r') is the ramping function defined as 

0 if <-T 

f{t') = < ( t' + T)/T - T < if < 0 (28) 

1 if > 0 

with T denoting the time constant, and < h(p / ) representing the velocity profile in cylindrical 
coordinates. Eq. (27) shows that the velocity distribution along the flight path (z-axis) 
is constant in the region of Li < z < L 2 . Also notice that the time constant T is 
approximately the rolling-up time of a vortex. The rolling-up distance for an elliptically 
loaded clean wing is at least between 2-3 aircraft spans downstream [28], thus on an 
average, 

T«2.5^ (29) 

where s is aircraft wing span, and U is the landing approach speed. 

With the definition of the ramping function f(t') in Eq. (28), the tensor, defined in 
Eq. (20), then becomes 

w,u) = ft4>'<S> 2 (p')[S(z' - Lj-Siz’ -L 2 )] 
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Eq. (30) shows that tensor A is an even function of the velocity profile, &(p'), therefore 
another vortex of the opposite sign will generate the same sound pressure, except an extra 
consideration of the distance to the observation point. Note that Eq. (27) represents a line 
vortex that is rolling up into a fully developed vortex in a longitudinal region bounded by 
Li < z < L 2 during time T, where L 2 — Li is determined by the coverage area of an array 
of acoustic wave sensors. Also note that the line vortex situated in a finite layer, so that 
the vortex held is localized and the radiation condition [29] is satisfied in the formulation. 
The assumption reflects the physical situation that local vortex sound is sampled by using 
a directional acoustic sensor, such as a phased microphone array as discussed in Sec. V-D. 
Upon introducing A in Eq. (30) into Eq. (24), the spectrum is found as 

P(f,u) = tvace{B(u)-!ft > —e : ’ ujr / Co J ft ft ■ ff < h 2 (p / ) 

• [S(z' - Li) - S(z' - L 2 )} e- j ^ f ' /Co df') (31) 

where 

B(u) = % [ — ^ fl - e~ jujT ) - 1 
K J T [ jTuj V > 

Furthermore, given the position of observation point f = — ry as shown in Fig. (1), and 
denoting f' = xp' coscp' + yp 1 sin (f/ + zz' and ft = —x sin d/ + ycoscj)', Eq. (31) is then 
obtained as 

P(f,u) = B{u)-^{L 2 -L x )ei“ r l c ° 

Cgirr 

• J cos 2 ^(//)e , '^ sh VW (32) 

Noticing that the integration with respect to the angle <fi gives 

J cos 2 ft e j ^ p smrf> dft = 7 r [J 0 (c op! j eft) + J 2 (up / c 0 )\ (33) 

where Jo and J 2 are the zeroth and second order Bessel functions of the Erst kind, respec- 
tively. In addition, by using the recursive formula [30] 

2jTi 

Jn-l(z) + Jn+ l{z) — —J n (z) (34) 

Eq. (32) is furthermore simplihed as 

°o 

P(f,u) = (L 2 ~ Xi) ft“ r ' Co f <$> 2 (p)Ji((jjp / cftdp 

c Q ur J 


( 35 ) 
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The above integration involves a single-fold integration in terms of radial distance p', 
therefore the computation is more efficient than the multiple-fold integration of the pre- 
vious expression. 


A. Example Vortex Profiles 


The two vortex profile models considered in this report are the Hallock-Burnham vor- 
tex [3] and the Lamb-Ossen vortex [3]. Both profile models are axisymmetric, with the 
former expressed as 


= 


P‘ 


12 


2np' p' 2 + p 2 

where T 0 is the circulation and p c is the vortex core radius, and the latter defined as 


(36) 


t(p') = 


2np' 


exp 


-1.2526 ( — 
\Pc , 


(37) 


The core radius, p c , defined as the distance between the locations of zero tangential 
velocity to where the maximum value occurs, is the same for both of the vortex profile 
models. However the peak values are different, e.g., 


for the Hallock-Burnham vortex and 


T 


o 


4vr p c 


JWi -L2526\ 

2tt Pc V > 

for the Lamb-Oseen vortex. For both of the vortex profile models selected, closed form 
solutions of the sound spectra are not achievable. The acoustic spectra of these vortex 
profile models are therefore obtained via numerical integration of Eq. (35). 


B. Validation Using Benchmark Vortex Profile 

For the purpose of validating the numerical integration program, a benchmark vortex 
velocity profile which allows a closed form solution to Eq. (35) is prescribed as 


<%') = 


47rp; 


2 pe Pc 


(38) 


with the core radius, p c , and the peak value, r o /(47rp c ), being the same as the Hallock- 
Burnham vortex, but with a faster decay for p > p c . With this definition of the vortex 
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profile, a closed form formula of the spectrum of sound pressure, with the derivation 
shown in detail in Appendix B, is obtained as 

P (39) 

where 


and 


A(u) = 


7T 


(a 2 + b 2 ) 
2an 

+■- — 

(\/a 2 + b 2 


(5«o + l) 


5a 2 7r 


4 (a 2 + 6 2 ) ; 
5a 3 7r 


4 (Va 2 + b 2 


; a D (5a 2 + 3) 

(a 2 + l) (5a 2 + 1 

(«o + l ) 2 


7 


a = 
b = 


2 C 0 
p c Ld 
L0jc 0 

2/Pc 


'1 + 


PcU 
2 Co 


(40) 


(41) 

(42) 

(43) 


It demonstrates that the vortex sound is independent on the sign of circulation, r o . Also 
notice that the spectrum of the vortex pressure is roughly proportional to T 2 and inversely 
proportional to p\ at high frequency. This implies that, taking into account that T a /p c is 
roughly a constant for current civil transport aircraft, a vortex with a large core (usually 
produced by larger aircraft), generates weaker sound in the high frequency regime. 

C. Numerical Results and Discussions 

The numerical evaluation of the spectrum expressed in Eq. (35) is then applied to a 
line vortex with the two commonly employed profile models indicated in Eqs. (36) and 
(37). As indicated earlier, a numerical code was developed to carry out the integration in 
Eq. (39), and the vortex profile of Eq. (38) was devised in order to provide a closed form 
solution for the purpose of validating the computer program. Characteristics of B747 
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TABLE I 

Constant Parameters for Various Aircraft 



Circulation 

Wing Span 

Landing Speed 


[m 2 /s] 

[m] 

[m/s] 

B747 

600 

60 

66 

B757 

360 

38 

66 


TABLE II 

Rolling-Up Time Constant and Core Radius 



Time Constant [s] 

Core Radius [nr] 

B747 

2.27 

4.67 

B757 

1.44 

2.88 


and B757 vortices generated during landing are used for this validation effort, whose 
circulation taken from [31] and core radius estimated from [32] 

Pc = 0.2^/hd (44) 

where s is aircraft wing span and U is the landing speed, are listed in Table I. From 
Eqs. (29) and (44), the rolling- up time constant T and core radius are estimated as shown 
in Table II. Other parameters adopted in the numerical calculation are air density [3], 
p 0 = 1.2 kg/m 3 , the sound speed, c a = 340 m/s, longitudinal window size L- 2 — L i = 20 m, 
and the distance between vortex and observation point, r = 60m(~ 200 ft). The choice 
of the particular distance represents the nominal aircraft altitude at the middle marker 
location, which has operational significance from wake turbulence perspectives. 

The numerical result of the pressure spectrum for the benchmark (BMK) vortex profile 
of Eq. (39) is compared against the closed form solution and shown in Fig. 3. Note that 
good agreement is found, where the Sound Pressure Level (SPL) in decibel scale, P r = 
20 log 10 (P/P D ), is computed using the Threshold of Hearing (i.e., P Q = 2 x 10 -5 N/m 2 ) as 
the reference (see Appendix C). The overall trend in the calculated spectra is similar to 
the recent phased microphone array measurements of [20] for rolling-up vortices generated 
during out-of-ground effect. It is further noted that the computed spectra do not suggest 
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the presence of distinct audible high frequency sound that has sometimes been reported 
in field observations. A candidate explanation is that since the prescribed vortex is of 
laminar state in a quiescent environment, mechanisms such as small scale turbulence and 
their interactions with the vortex which could manifest themselves in higher frequency 
sound, were not considered. Otherwise stated, the computed spectra is based on the 
aspect of vortex dynamics that is always present and robust, and can be considered as 
a baseline spectra. The higher frequency sound sometimes heard when vortices are near 
the ground is speculated to be caused by the turbulence generated from vortex-ground 
interactions. This would explain the unpredictable nature in the audible portion of the 
vortex sound [4], 

The good agreement provided the necessary confidence to employ the numerical code 
with the more realistic vortex models of Eqs. (36) and (37). Characteristics of two specific 
aircraft types, B747 and B757, are listed in Table I. ft is worth noting again from Eq. (39) 
that the core radius is an important parameter. However, a satisfactory universal scaling 
scheme for the core size relative to aircraft and atmospheric parameters is still an area 
of research. Core radii listed in Table I may be over-estimated, but can still be taken as 
values reflecting the relative sizes of the two example aircraft. 

Velocity profile models for the Hallock-Burnham, Lamb-Oseen and the benchmark pro- 
file (38) leading to a closed solution of Eq. (39), generated by B747 and B757, are il- 
lustrated in Fig. 2. Note that, in all of these profile models, the peak velocities appear 
at p = p c . In addition, the peak values are the same in comparison between B747 and 
B757, owning to the fact that the peak value is proportional to the ratio of T 0 /p c which 
is constant in the case presented. Also notice that the benchmark profile has the same 
peak value as the Hallock-Burnham vortex, but decays more rapidly after p > p c than 
both Hallock-Burnham and Lamb-Oseen vortices, ft also shows that the vortex profile for 
larger aircraft (e.g., B747) spreads out wider than that for smaller aircraft (e.g., B757). 

Fig. (3) also shows the numerical results of sound pressures for Hallock-Burnham and 
Lamb-Oseen vortices generated by B747 and B757, respectively, at the distance of 60 m 
(~ 200 ft). The integration in Eq. (35) is performed by using trapezoidal method. In 
comparison with vortices with different core sizes, it is found that the vortex with faster 
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Fig. 2. Velocity profiles for the benchmark (BMK), Hallock-Burnham, and Lamb-Oseen vortices gen- 
erated by B747 (upper plot, r o = 600m 2 /s, p c = 4.67m) and B757 (lower plot, r o = 360m 2 /s, 
p c = 2.88m), respectively. 

decaying velocity produces a wider frequency band for its sound pressure. For example, 
the B757 generates a wider spectrum of vortex sound than the B747 does. On the other 
hand, as expected, larger aircraft (B747) create higher vortex sound around the peak 
value. 

D. Gain Pattern of Sensor Array 

The factor L 2 — L\ in Eqs. (35) and (39) indicates the vortex length from which the 
sound is detected. Physically, it can be implemented by using a linear phased microphone 
array deployed under the vortex as illustrated in Fig. 4. 
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Fig. 3. Spectra of sound pressure for Hallock-Burnham (HB), Lamb-Oseen (LO), and benchmark (BMK) 
vortices, with constant velocity distribution along flight path, generated by B747 and B757, respec- 
tively, at distance r = 60m(~ 200 ft). For B747, the circulation is r o = 600m 2 /s, and the core 
radius is p c = 4.67 m; for B757, the circulation is r o = 360m 2 /s, and the core radius is p c = 2.88 m. 

Denoting the locations of the vortex source and the n-th sensor as z and rid, respectively, 
the distance between the vortex and the sensor can be expressed as 

L n (z) = yj \{z - nd) 2 + r 2 ) (45) 

In addition, by assuming that the total number of sensors in the array is N, the spacing 
between sensors is d, and the height of the vortex region is r, then the output, which is 
the coherent sum of the signals from all sensors, is obtained as 

n =(iv— 1)/2 

S(f,z)= J2 P(f) ex P (jkL n (z)) 

n =-(N- 1)/2 


( 46 ) 
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Fig. 4. Layout of acoustic sensor array along aircraft flight path. 

where k is the wavenumber, k = 2nfjc 0 , with / and c Q denoting the frequency and speed 
of sound in free space, respectively, and P(f) is the Fourier transform of the sound signal 
detected by a single microphone sensor. Finally the normalized gain pattern of the sensor 
array is calculated by 

G(M = 2oi ° g (!M) (47) 

Fig. (5) shows the numerical result of the gain pattern for the array operated at fre- 
quency / = 50 Hz, where the spacing between elements is d — 1 m, the number of elements 
N = 19, and the height of the vortex region is r = 60 m. Note that the region for the 
gain higher than -3 dB is —10 m < z < 10 m, implying that only the vortex in the range 
of [—10 m, 10 m] contributes significantly to the output after the combination of signals 
from the sensor array, as assumed in the numerical simulations for the vortex sound spec- 
trum performed in this report. Another usefulness of using the sensor array is to enhance 
the power of the sound spectrum. If M arrays are used with each array possessing N 
elements, the output power of the spectrum will gain approximately 20 log(MTV) dB. For 
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Fig. 5. Gain pattern of acoustic sensor array with TV = 19 elements along aircraft flight path for vortex 
sound frequency of / = 50 Hz. The array spacing is d = 1 m. 

example, if M — 5 and N = 19, then the gain will be 40 dB. 

VI. Summary and Closing Remarks 

The spectral characteristics of sound generated by the rolling-up process of wake vortices 
at out-of-ground effect altitude have been analytically examined. The acoustic signature 
was derived from a tensor analysis on the acoustic analogy equation, which ultimately 
lead to a general analytical expression of the spectra. The influence of the core radius, 
vortex velocity profile, and circulation are also examined by performing calculations for 
two representative aircraft, namely, B747 and B757. The vortex sound spectrum was 
found to be of broadband nature, with a characteristically higher sound level located in 
region below 100 Hz. The SPL drops off very rapidly after the aforementioned region 
with increasing frequency. For example, the calculations revealed that at an altitude 
of 60 m (~ 200 ft) away from the vortex, the sound pressure is as high as 60 dB with 
respect to the Threshold of Hearing, which is equivalent to listening to softly spoken 
human voice from a distance of 1 ft (0.3m). While for frequency higher than 100 Hz, 
the sound pressure is lower than the TOH, indicating that it would not be audible in a 
typically well trafficked airport environment. It is important to note that turbulent flow 
mechanisms around the cores have not been modeled in this analysis. The computed 
spectra is therefore representative of the baseline characteristic from a line vortex during 
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rolling-up. The high energy content situated around the infrasound regime may render 
itself conductive to wake monitoring, as the atmospheric attenuation effect is less severe. 
However, it may be more practical from a sensor perspective to monitor frequencies higher 
than the infrasound regime. Finally, as the speed of sound is relatively slow (compared to 
light) in reaching a detector which translates to a level of uncertainty in the detected wake 
locations, the passive acoustic scheme of wake monitoring, if can be proven to perform 
reliably, is envisioned to be a sensor with a practical range less than a mile. 


VII. Appendices 

A. Derivation of Eq. (22) 

The function F(r', u) in Eq. (19) is the Fourier transform of /(f',0), written as 

F(f',u) = J f(r',t')e ju}t 'dt' 

= V ■ V ■ J [u(f', t')u(f\ t')] e juJt 'dt' 


By using the identity 


0 jui\r-r’\lc 0 


(f)V ■ a = V ; • (0a) — a ■ V0 


(48) 


(49) 


where 0 = — — — — and a = V 7 • A(r',u) with the tensor A expressed explicitly in 
4-7T |r — r'\ 

Eq. (20), and ignoring the surface integral term assuming the vortex disappears at infinite 
distance [29] 1 , the spectrum is found as 


P(r,w) = - Po j (V ■ .4(f',w)) ■ V<pdf' 


(50) 


Furthermore, applying the identity 

(V • I) • b = V ■ (A ■ b) - trace (a ■ Vbj (51) 

where b = V ; 0, again ignoring the surface integral and notifying the symmetry property 
of the tensor, A = A , the compact form of the spectrum is then obtained as shown in 
Eq. (22). 


lr The Divergence Theorm is applicable both for vector and tensor A with the form J V ■ Adr = ^ ds-A. 
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B. Closed Form Spectrum for Benchmark Vortex 

Introducing the velocity profile of the benchmark vortex, Eq. (38), into the integral form 
of spectrum, Eq. (32), after conducting the integration over p', the following expression 
is obtained, 


P(r, u) = j 


. 3 p Q e 2 ( T 0 \ 


2 c IpW 


2ir 


(^) ( L 2 -U)e j “ r/C ° 

cos 2 0' 


J (uj/ c Q sin 0' + j2/p c ) 4 

Note that the integral term has the general form 


d(j)' 


2ir 


COS 2 0' 


A (c j) = 


rdf)' 


(52) 


(53) 


4 1 

^ (asin0' + j5) 

where a = c o/c Q and b = 2/ p c . By making change of the variable, z = e^' , Eq. (53) 
becomes a contour integral along a unit circle on the complex plane, 


z{z 2 + lf 


A(u) = — j . 

ja J + (, < 

where the integrand can be expanded as 


b Va 2 +b 2 


rdz 


(54) 


Integrand = 
-4 


^ 2 + iy 


1 1 

u + 


(a Q - fd 0 ) { (z - a Q ) (z - (3 0 ) 

1 1 4 1 1 

(a Q - Po ) (z - a 0 ) 3 (a Q - p o ) (z - p o ) 3 

11 11 

+ 10 o o + 10 o o 

(ck 0 - Po) (z - a 0 ) (a Q - Po) (z - Pop 

-20 1 oT 1 r +20 1 77 1 — l 

(a Q - Po) 3 (z ~ a Q ) (a„ - P 0 ) 3 (z ~ Po) J 


with a Q and P a expressed, respectively, as, 


„ b , 

«o, Po — - zz 

a a 


Then using the Cauchy’s integral formula 




n\ r f(z) 


-dz 


(55) 


(56) 


(57) 


2iJ J (z-a 0 ) n+1 

the integration of Eq. (54) is found in a closed form as Eq. (40). Finally the closed form 
solution of the sound spectrum for the benchmark vortex profile is obtained as Eq. (39). 



23 


C. Vortex Sound Pressure Level vs. Threshold of Hearing 


The vortex sound, P, is expressed as the Sound Pressure Level (SPL) in decibel scale 
with respect to the pressure of the Threshold of Hearing (TOH). The intensity of the 
TOH is I Q = 10 _12 W/m 2 for a single tone at 1 kHz [33], thus the pressure of the TOH is 
P Q = 2.0 x 10 _5 N/m 2 according to the relation [25] 

pi_ 

PoCo 


L = 


(58) 


A 0-dB SPL represents the level of sound pressure at P = P Q , and a 74-dB SPL is 
equivalent to the sound pressure of human normal speech at 1-foot (0.3 m) distance. 
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